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Abstract 

In terms of a simple holographic model, we study the absorption cross section 
and the shear viscosity of a pure Yang- Mills field at low temperature where the 
system is in the confinement phase. Then we expect that the glueball states are 
the dominant modes in this phase. In our holographic model an infrared cutoff 
r m is introduced as a parameter which fixes the lowest mass of the glueball. As a 
result the critical temperature of gluon confinement is estimated to be T c ~ 127 
MeV. For T < T c , we find that both the absorption cross section and the shear 
viscosity are independent of the temperature. Their values are frozen at the 
values corresponding to the critical point, for < T < T c . We discuss this 
behavior by considering the glueball mass and its temperature dependence. 
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1 Introduction 



Recently, Policastro, Son and Starinets pQ have given a nice application of holographic 
calculations for the shear viscosity. This viscosity could be observed in the hot QGP 
phase of QCD. The calculation is based on the gravity/gauge correspondence [2l[3jlU[5], 
and assumes a strongly coupled 4d YM theory. 

With regard to experiments the situation is however unfortunate, since the quark- 
gluon plasma one hopes to create in heavy-ion experiments has relatively low temper- 
ature at which the coupling constant is strong and the perturbation theory in QCD 
works poorly. So we expect that an approach in terms of holographic models, utilizing 
the gravity/gauge correspondence, would be very useful for the understanding of the 
various experimental result in this region. 

Recently, we have found that the holographic model available at high temperature 
can be extended also to the low temperature quark confinement phase by introducing 
an infrared cutoff [HllTj. This approximate model could give successful results for the 
thermal properties, masses of light mesons, decay constants etc, in the low temperature 
confinement phase of 4d YM theory. In the confinement phase, the physical states of 
QCD would be hadrons. Then we would be able to see the thermodynamical properties 
of the hadron gas by applying our approximate model to the 4d YM theory in the low 
temperature phase. 

In the present paper we study the shear viscosity, which is related to the absorption 
cross section of gravitons []]. We first study the absorption cross section of gravitons 
falling perpendicularly onto D3-branes in the confinement phase. The calculation is 
performed in a way similar to that of Refs. [EJ [9]. Then we extend the calculation to 
the finite temperature case in order to see the temperature dependence of the shear 
viscosity of hadron gas, i.e. the glueball gas in the present case. 



2 Absorption by deformed D3 Brane 



The shear viscosity is related to the absorption cross section, so we firstly study this 
quantity at zero temperature to understand the situation in the confinement phase. 
The D3-brane metric can be written as 



ds 2 = H- 1 ' 2 (-dt 2 + dx\ + dx\ + dxf) + H 1 ' 2 (dr 2 + r 2 dnj) , 



where 



H 



The s-wave of a minimally coupled massless scalar satisfies 

(loR) 4 ' 
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where p = tor. For small uR the absorption cross section was obtained by Klebanov 
[8] by solving the matching problem of an approximate solution in the inner region and 
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the outer region. The following result was obtained, 



0"3-brane — R 



(2) 



And this is also precisely given from 4d A/"=4 SYM theory as a three point vertex of 
the bulk field and the fields of SYM theory which are given by the Born-Infeld action. 

In this case, the 4d J\f—4 SYM theory is described by string theory or by supergrav- 
ity near the horizon of the D3-branes. In this limit, the geometry is approximated by 
AdS*5 x S 5 and the equation of motion for the dilaton (0) or for the graviton polarized 
parallel to the D3 brane is written as 



where M g denotes the 4d mass of the dilaton. However, for AdS^ x S 5 , there is no 
finite M g which gives a normalizable wave function since J\f=4 SYM is not a confining 
theory. To obtain a confining YM theory, the bulk geometry should be deformed to 
an appropriate form. In this case, however, the geometry has a complicated form 
especially in the infrared region. 

A simple and tractable approximation for such a confining geometry can be given 
by introducing an infrared cutoff at an appropriate point of the coordinate r, say 
r = r m . In this case, the geometry AdS§ x S* 5 is used in the region r m < r < oo. This 
approximation would be justified since the main deformation of the geometry is in the 
infrared region, r < r m , and this region is cut off now. In fact, we can obtain a discrete 
mass spectrum for the glueball by solving (j3J) and imposing the boundary condition 
dr<fi\r m = 0. When we adopt the lattice result for the glueball mass, M g ~ 1.5 GeV, 
the IR cutoff is set at r m ~ 0.4 in units where R = 1. Here we used the data in [TO] 
since the temperature dependence of the glueball mass is given in this reference. We 
see that it is consistent with our analysis given in Fig{TJ 

Now we turn to the absorption cross section for the deformed D3 branes, whose 
horizon limit background is given by the AdS§ x S 5 with IR cutoff. Then the back- 
ground, which describes also off the horizon limit, can be written as ([2]) with the IR 
cutoff r m . So the matching problem is solved in a similar way as that developed by 
Klebanov. The only difference is the existence of the IR cutoff, and the equation is 
restricted to the region r m < r < oo. We call the region near r m the inner region. 

While the wave equation is solved in terms of Bessel functions at small p [8], we 
here solve the equation at some small but finite r, namely at the IR cutoff r m (< R). 
To study the wave function in this region, Eq. ([1]) is solved by rewriting it according to 
as 




(3) 



d 2 T + u 2 m 0(r) = O, r 
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(4) 
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The solution of this form of equation can be matched to the outer solution. Then, near 
r m , the incoming wave is given as 

0(r) ^ = e-*+» T (5) 
and this is approximated for small r as 

4>(t) in = 1 - iu m T (6) 

In the outer region, (DO) is estimated by substituting <p = p~ 5/ ' 2 V ; - Then we get 

[^-v + 1 + V- w = - (7) 

Now the last term is negligible for p ^> (ujR) 2 , where (J7J) is solved in terms of Bessel 
functions as 

(Pom = Ap- 2 [J 2 (p)+btN 2 (p)} , (8) 
J and iV being the Bessel and Neumann functions. 

From the matching of the solutions, (jHJ) and (jHJ), in the overlapping region of r, we 
obtain A = 1 and 6 = 0. Then 
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'out 



p-V 2 (p) • (9) 



In terms of the formula for the invariant flux F, 

F = ±:V=gg rr ((l>*d r cp-c.c.), (10) 

the fluxes of the incoming wave are obtained as 

-dinner = ^9m, (H) 

g m = r 5 m ^/l + (R/r m y 

near r m , and 

^outer = — & 4 - (12) 

7T 

for the outer region. Then the absorption probability at r m is given as 

I -dinner i ^ 5 / 1 o\ 

" outer 

We notice that the power of u is different from the case of D3 background without the 
IR cutoff. 

The absorption cross-section a is related to the s-wave absorption probability by 
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where 

27r -2±i 
Q n = — 7 \~ 

r(s±i) 

is the volume of a unit D-dimensional sphere. Thus, for the 3-brane we find the 
following cross section at the infrared cutoff r m , 

cm = TT 3 g m • (14) 

We notice that this result is independent of to contrary to the case of the background 
without the cutoff, and it is proportional to the area of the five dimensional sphere 
with radius r m . So the situation is similar to the absorption by a black body in the six 
dimensional space. 

In the case of background without the cutoff, the cross section is explained by the 
decay amplitude of the incoming scalar to two gluons or other fields included in the 
SYM theory [HI [TTJ [T2]. In this case, the semi classical calculation is useful because 
of the conformal invariance of the YM theory. And the gluons are not bounded to 
glueballs in spite of the interactions with a strong coupling constant. In the present 
case, however, the conformal invariance of the gauge theory is of course lost and the 
gluons are confined and bounded to glueball states. Then we expect that the incoming 
scalar of energy u, which is in the range M g < uj < 2M g , will form a glueball after it 
coupled to the gluon pair. For uj > 2M g , there appears the possibility of pair production 
of glueballs, and it would be possible to calculate the pair production cross section as 
in the YM case. But the interaction between the glueballs is complicated and strong. 
It would be difficult to get some useful results from the 4d field theory side. This point 
will be discussed more in the future. 



3 Extension to finite temperature 



For the finite temperature gauge theory, the corresponding gravity background solution 
is given by the nonextremal black three-brane. It has the form [131 EE 



ds 2 = H-V 2 (r)[-f{r)d1? + dx 2 ] 

+H 1/2 (r)[f- 1 (r)dr 2 + r 2 dnl] J (15) 

where as before H(r) = 1 + i? 4 /r 4 and f(r) = 1 — r 4 /r 4 . The Hawking temperature of 
this metric is 

T= J j j5 , (16) 

and the metric of the previous section is obtained for T = 0. It is well known that 
gravity under this background in the near-horizon limit of D3 branes describes the 
high temperature deconfmement phase of the YM theory. On the other hand, we have 
proposed that the low temperature confinement phase could be obtained by introducing 



*By absorption cross section we will consistently mean the cross section per unit longitudinal 
volume of the brane. 
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the IR cutoff, r m , in this holographic model as r < r m (< R) P [7]. When the two 
phases are separated at r = r m , the critical temperature T c is given as 

T c = ~ 127 MeV . (17) 

Next we consider the s-wave radial equation for a minimally coupled scalar (such 
as the graviton polarized parallel to the brane), d^^—gg^d^cj)) = 0. For the metric 
(fT5l) this equation acquires the form 

d P { P *fd p <P) + ^0 = . (18) 

We solve this equation in two phases, one at low and one at high temperature, and 
give the absorption cross sections in each phase. 



3.1 Low temperature phase T < T c 

By introducing the IR cutoff r < r m (< R), the low temperature phase is obtained 
for the range of r m < r < oo in terms of the metric ( 1T51) . Here, we consider r m to be 
independent of the temperature or of tq. We estimate the absorption cross section at 
the IR cutoff as in the previous section. Near r = r m , Eq. (fT8|) is rewritten as 

a!+i4,]0(f) = o, (w) 

9f 1 (20) 



dr p 5 f ' 

where we notice u> m = p b m \Jl + (uRy/p^, p m = ur m . Then at the IR cutoff, the 
incoming wave is given as 

0(r) in = e~^ f , (21) 
and this is approximated for small f as 

1 - iu m f. (22) 



m 



On the other hand, in the region of large p, we can approximate H ~ 1 and / 
then we arrive at the same form of solution as in the previous section, 

<j) = Ap- 2 [j 2 (p) +biN 2 (p) 

As in the previous section, the outer solution is given as 

0out = ^P~ 2 J 2 (p) > 

and this matches with (|2T|) at an appropriate region near r m . Then we arrive at the 
same result for the fluxes of the incoming wave as that obtained for T = 0, 

-dinner = ^Qmi (23) 
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where g m = r^y 1 + {R/r m ) i as shown above. We should notice here that the factor 
f(p) has disappeared in the final result of (I23I) due to the exact cancellation between 
the wave function and the measure in (fTOj) . As a result, F inner is independent of the 
temperature. 

Now, -F outcr is the same as given in (fl2|) at zero temperature. Then all other quan- 
tities are the same as in the zero temperature case. The absorption probability at r m 
is given as 

,-P, I -firmer i ^ 5 fnA\ 

V = \p \ = ^9m, (24) 

-^outer "J^ 

and we obtain the absorption cross section at low temperature, 

a 3°^brane = ^ 9m ■ (25) 

This result is seen to be independent of the temperature and it has the same form as 
that obtained at zero temperature. This point can be understood as follows. Since 
T < T c << M g , the thermal fluctuations of glueballs would be suppressed because we 
need large energy, at least uo > M g , to obtain a definite thermal effect, which should in 
principle be expressed by some kind of temperature dependence, according to a thermal 
field theory for the glueball gas. 

In the above discussion we neglected the temperature dependence of the glueball 
mass M g (T), but it will be important. It is estimated here by solving the equation of 
motion of the dilaton near the D3 horizon limit of (jT5l) . 




Fig. 1: Glueball mass vs temperature, for r m = 0.4 and R — 1. 



The equation to be solved is 

(d 2 r + d r (log(r 5 /(r))) d r + M g 2 l + ( ^ /r)4 ) = 0, (26) 
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and the mass eigenvalue M g is obtained by finding the normalizable solutions of this 
equation. The numerical results are given in FigJT]for r m = 0.4 and R = 1. We observe 
its typical T-dependence, M g — > for T — > T c , and this behavior can be understood by 
the last term in the equation. The similar behaviour has been seen also in [BJ [T] and 
it is consistent with the lattice analysis [TU] as pointed out above. From these results 
we can say that the above T-independence of (1251) would be modified near T c since M g 
becomes small at this point and the probability of pair production of glueballs can not 
be neglected in this region. 

It is obvious that the temperature dependence of the mass comes from the factor 
/ in the equation. On the other hand, the factor / is absent in ( 1251) . Therefore, the 
temperature dependence has to be included in the parameter r m itself, when o"3_brane is 
temperature dependent. In other words, r m will vary with temperature, at least near 
T = T c . This point will be discussed in the next sub-section. 



3.2 High temperature phase 

At high temperature, r m < ro, we consider the region of ro < r < oo where we can see 
the deconfinement phase of the gauge theory. Here we estimate the absorption at the 
horizon tq. Near r = ro, Eq. ( Jl8|) has the form 



'," + -i— + 



f^V^ = 0. (27) 



/• — /•(i \4nT J (r — ro) 2 

The matching problem has been solved as in pQ, and we obtain the absorption proba- 
bility 

,-!-, I -^near horizon i ^53 r>2 

V = = — l) rZR 

± outer 

and the cross section 

^brlne = ^fgi? . (28) 

Then the cross section increases like T 3 in the high temperature QGP phase. 



On the other hand, the low temperature result Eq. fl25|) is rewritten for the case 
R/r m » 1 as 

4° W bTane ~ ^R 2 , (29) 

and we see that this is equivalent to fl28l) for ro = r m . This implies that the thermal 
activity is frozen by the value at T = T c for T < T c since the glueball mass is much 
larger than T c . All the degrees of freedom contribute to the thermal property for T > T c 
and we obtain the temperature dependence of f[2"B"j) . 

As for the IR cutoff parameter r m , it has been introduced by hand as a constant 
which is independent of ro- However, it may have a temperature dependence. In such 
a case, we can refine the description and write r m = r m (ro), imposing the following 
condition ^ 

r m {T c ) = r c , T c = —j^. (30) 
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Then the value of T c is also changed to a larger value than the 127 MeV given above. 
In fact, it is estimated as T c = 270 MeV in the lattice gauge theory [10]. When we take 
this larger value, we obtain r c ~ 2r m (T = 0). And we thus obtain the temperature 
dependence of 0" 3 

*KLeCTe) = 8^_ W b Tane(0) • (31) 

However, this analysis should be studied in more detail and it is outside the scope of 
the present work. 



4 Shear viscosity 

We briefly review the notion of viscosity in the context of finite-temperature field 
theory. Consider a plasma slightly out of equilibrium, so that there is local thermal 
equilibrium everywhere but the temperature and the mean velocity slowly vary in space. 
We define, at any point, the local rest frame as the one where the three-momentum 
density vanishes: T 0i = 0. The stress tensor, in this frame, is given by the constitutive 
relation, 

Tij = 5 i:j p - 77 (diUj + djUi - ^SijdkUk 

-(SijdkUk , (32) 

where Ui is the flow velocity, p is the pressure, and rj and ( are, by definition, the shear 
and bulk viscosities respectively. In conformal field theories like the jV=4 SYM theory, 
the energy momentum tensor is traceless, T 7 ^ = 0, so e = T 00 = 3p and the bulk 
viscosity vanishes identically, ( = 0. 

All kinetic coefficients can be expressed, through Kubo relations, as the correlation 
functions of the corresponding currents [15]. For the shear viscosity, the correlator is 
that of the stress tensor, 

■niu) = T^Jdtd* ^ < [T*y(t> x),T xy (0, 0)] > 

= ^[GaM-GrH, (33) 

where the average < . . . > is taken in the equilibrium thermal ensemble, and and 
Gr are the advanced and retarded Green functions of T xy , respectively. In Eq. f l3"3"j) . 
the Green functions are computed at zero spatial momentum. Though Eq. (1331) can, in 
principle, be used to compute the viscosity in weakly coupled field theories, this direct 
method is usually very cumbersome, since it requires resummation of an infinite series 
of Feynman graphs. This calculation has been explicitly carried out only for scalar 
theories [17J. A more practical method is to use the kinetic Boltzmann equation, 
which gives the same results as the diagrammatic approach [16J. 

The static shear viscosity is defined in the limit of uj —>■ of the correlator of energy 
momentum tensor, 

77(0) = \imj- f dtdxe iuJt ([T xy (t,x),T xy (0,0)}) 



1 

2~K 



2 *(0) (34) 
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In the high temperature phase, the viscosity is given from (1281) as [lj 



W = — T 3 . (35) 



On the other hand, in the low temperature phase, we obtain by using Eqs. fl25|) and 
dMD, 



ttN 2 ( r m \ 3 (r, x 2 



m 



TtR 2 ) \R) \ 



1 (-)'• (36) 



When we do not consider z m as the T-dependent parameter, this low temperature 
viscosity is independent of the temperature. When we restrict ourselves to the case 
r m /R « 1, it is approximated as 

TtN 2 ( r m Y ttN 2 3 
5- -7^ = -^^ c • (37) 



8 \nR 2 

It coinsides with f l35l) at T = T c . 

This result seems to be reasonable. Since the lowest state of the glueball has a 
large mass (~ 1.5 GeV) and its gas would not be thermodynamically active below 
T c ~ 127 MeV, we expect its T-dependence to be negligible. When we adopt the 
critical temperature as T c ~ 270 MeV [TU] , however, the cross section changes as given 
in (!3T|) . As a result, we expect the T-dependence of the viscosity for YM theory to be 
as shown in FigfJl 



?7/?7c 
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Fig. 2: 77/^ is shown, where r] c = Z ^T C 3 . Temperature is scaled by T c . The dashed 
line shows the case of T-dependent r m with T c = 270 MeV. 
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5 Summary 



In this paper we have studied the absorption cross section and viscosity of pure Yang- 
Mills theory at low temperature where the theory is in the confinement phase. The 
gluons are confined and we expect that the glueball states are the dominant modes 
in this phase. In our holographic model, the infrared cutoff r m is introduced as a 
parameter, and it is fixed by lowest mass of the glueball which is estimated to be 
about M g = 1.5 GeV. Using this r m , the critical temperature of gluon confinement is 
estimated to be T c ~ 127 MeV. 

We expect that the thermal properties of the YM theory below T c will be determined 
by the glueball gas. However, for T < T c the glueball mass is much larger than the 
temperature, then the glueball gas will not show any active temperature dependence. 
As expected, we find that both the absorption cross section and the shear viscosity 
are independent of the temperature. The value of the shear viscosity for T < T c 
is frozen at the value corresponding to the critical point up to T = 0. The only 
possible temperature dependence of the viscosity will be seen through the temperature 
dependence of the cutoff itself as r m = r m {T). This is an interesting point, to which 
we will return in the future. But it is an open problem at the present stage. 

We have not considered the light flavor quarks and the light mesons constructed by 
them. If we take into account those mesons, the thermal properties at low temperature 
will be changed due to the light flavored mesons [18]. This point is also an open 
problem at the moment. 

Finally, let us make some remarks from a fundamental viewpoint on possible phys- 
ical interpretations of the quantity r m . From the above formalism it follows that r m 
can be looked upon in two different ways. First, r m can be considered to play the role 
of an extra horizon, in addition to the conventional horizon r$ that is present at finite 
temperature. While Policastro et al. [1] got cr{yj) ~ ui 3 for small u at T = and thus 
cr(0) = 0, we got <tir to be finite at the infrared cutoff (Eq. (THl) ). at zero temperature. 
This is to be compared with the result cr(0) = tt 3 TqR 2 in [TJ, which is finite at finite 
temperature. Secondly, our result for the shear viscosity at high and low temperature 
is given by Eqs. (l3"5"j) and (13T)j) . Here, the behaviour of rj L is seen to be influenced by r m . 
It is of interest to compare also this result with that obtained in pQ: rj = f(g 2 N)N 2 T 3 , 
where f(x) is a function satisfying f(x) ~ x~ 2 ln~ 1 (l/x) when i«l and f(x) = tt/8 
when x ^> 1. We thus see that the quantity r m can manifest itself in different ways 
physically. It hardly seems possible to associate r m simply with a unique conventional 
physical quantity. 
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